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' A survey is presented of various aspects of the interaction of charged 

"c^ ' particles with sohds. In the framework of many-body perturbation 

theory, we study the nonhnear interaction of charged particles with a 
free gas of interacting electrons; in particular, nonlinear corrections to 
^ ' the stopping power are analyzed, and special emphasis is made on the 

O ■ separate contributions that are originated in the excitation of either 

electron-hole pairs, single plasmons, or double plasmons. Ab initio 
calculations of the electronic energy loss of ions moving in real solids 
^ ' are also presented, and the energy loss of charged particles interacting 

^ . with simple metal surfaces is addressed. 
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1 Introduction 



A quantitative description of the interaction of charged particles with sohds 
is of basic importance in many different theoretical and applied areas pp. 
When a moving charged particle penetrates a solid material, it may lose 
energy to the medium through various elastic and inelastic collision processes 
that are based on electron excitation and nuclear recoil motion in the solid. 
While energy losses due to nuclear recoil may become dominant at very low 
energies of the projectile j2], in the case of electrons or ions moving with 
nonrelativistic velocities that are comparable to the mean speed of electrons 
in the solid the most significant energy losses are due to the generation of 
electronic excitations, such as electron-hole (e-h) pairs, collective excitations, 
i.e., plasmons, and inner-shell excitations and ionizations. 

In 1913 Niels Bohr published a seminal paper on the energy loss of charged 
particles interacting with electrons bound to atoms ^3 , which laid the ground 
for Bethe's quantum theory of stopping power jlj. The many-body prob- 
lem of interacting conduction electrons in metals was investigated by Bohm 
and Pines [5^ in the so-called jellium model, by simply replacing the ionic 
lattice of the solid by a homogeneous background which serves to provide 
neutrality to the system. Within linear-response theory, the electronic re- 
sponse of conduction electrons to external charged particles is determined 
by a wavevector- and frequency-dependent longitudinal dielectric function. 
In the self-consistent-field, or random-phase, approximation (RPA), the di- 
electric function of interacting free electrons was derived by Lindhard [H], 
and, subsequently, a number of workers have derived alternative dielectric 
functions that incorporate various many-body higher-order local-field correc- 
tions [3 IB] and band-structure effects p. The effect of dissipative processes 
occurring in a real metal and conversion of plasmons into multiple e-h pairs 
has been allowed phenomenologically by including a damping coefficient in 
the dielectric function [TUj. 

The validity of linear-response theory, which treats the perturbing poten- 
tial to lowest order, is not obvious a priori. While lowest-order perturbation 
theory leads to energy losses that are proportional to the square of the pro- 
jectile charge Zie, the energy loss of either positive and negative pions [TT] 
or protons and antiprotons ^21 ^| is known to exhibit a measurable de- 
pendence on the sign of the charge ^1 1^ HHI- Experimentally observed 
nonlinear double-plasmon excitations ^3 ^] are also beyond the realm of 
standard linear-response theory fOJ 1201 1^ , nonlinearities may play an im- 
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portant role in the electronic wake generated by moving ions in solids [^ 1^. 
and lowest-order perturbation theory breaks down when the projectile is ca- 
pable of carrying bound electrons with it P5!- 

The first full nonlinear calculation of the electronic stopping power of an 
electron gas was performed by Echenique et al. in the low- velocity limit. 
They used a scattering approach to the stopping power, and the scattering 
cross sections were calculated for a statically screened potential which was 
determined self-consistently by using density-functional theory (DFT) |2E]- 
These static-screening calculations were then extended to projectile velocities 
approaching the Fermi velocity of the target j2Zj- Alternatively, an effective 
charge can be assigned to the projectile [2H] and nonlinearities can then be 
investigated within quadratic-response theory, thereby providing results for 
arbitrary projectile velocities [221 • 

To calculate the electronic stopping power of real solids from the knowl- 
edge of the full band structure of the solid and the corresponding Bloch 
eigenfunctions and eigenvalues is a laborious problem. Hence, early theoret- 
ical investigations were based on semiempirical treatments of the electronic 
excitations in the sohd [^01 IHH IS21 ISHl Ell- Attemps to introduce the full 
electronic band structure in the evaluation of the electronic stopping power 
of alkaline metals for low projectile velocities include a one-band calcula- 
tion as well as a calculation based on a linear combination of atomic 
orbitals The low- velocity limit was also investigated, in the case of sili- 
con, on the basis of a static treatment of the density response of the solid [371 • 
Ab initio band-structure calculations that are based on a full evaluation of 
the dynamical density-response of the solid have been carried out only very 
recently gHl 123110]. 

In this paper, we summarize recent investigations on the impact of nonlin- 
ear, band-structure, and surface effects in the interaction of charged particles 
with solids. We present general procedures to calculate, within many-body 
perturbation theory, the nonlinear potential induced by moving ions in an 
interacting free-electron gas (PEG), Z\ contributions to the stopping power 
of a PEG, and double-plasmon excitation probabilities. Self-consistent cal- 
culations of the energy-loss spectra of charged particles moving near a metal 
surface are also presented. Pinally, we consider the electronic stopping power 
of valence electrons in real solids, which we evaluate within random condi- 
tions. 

Unless otherwise is stated, atomic units are used throughout this paper, 
i.e., = h = rrtf. = 1. The atomic unit of length is the Bohr radius. 
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o-o = h'^/fnl = 0.529 A, the atomic unit of energy is the Hartree, 1 Hartree = 
e^/ao = 27.2 eV, and the atomic unit of velocity is the Bohr velocity, Vq = 
ac = 2.19 X lO^cms"^, a and c being the fine structure constant and the 
velocity of light, respectively. 



2 Theory 

We consider a recoiless particle of charge Zi ^T] moving in an arbitrary 
inhomogeneous electron system with nonrelativistic velocity v, for which 
retardation effects and radiation losses can be neglected. The energy that 
the probe particle loses per unit time due to electronic excitations in the 
medium can be written as 142 



dE r ,,,, .dV'^\v,t) 



(1) 



where p^^*(r, t) represents the probe-particle charge density 

p^^'*(r,t) = Zi5(r-ro-vt) (2) 

and V^*"-'^(r, t) is the induced potential. Keeping terms of first and second or- 
der in the external perturbation, time-dependent perturbation theory yields 

r-(r,t) = fdr' Tdt' T'^e-^^-^UdrJdr, 

J J — OO J — CO 2?? J J 

^ dtu' 



X 



X 



X{ri,r3;uj)+ dr" dt" 

J J —CO 

I dr2 Y{r^, rs, r^; u, cu') t;(r2, r") p^^*(r", t" 



v[r, ri) 

i{uj+ij'){t'-t") 



(3) 



where f(r, r') = l/|r — r'| is the bare Coulomb interaction, and 

and Y{yi, r2, ra; tu, uj') are the so-called linear and quadratic density-response 

functions of the electron system: 



and 



F(ri,r2,r3;a;,a;') 



(ri)]o. 



nO 



[r2)]o„ [P(ri)] 



raO 



U -UJnO + iV + UJnO + if] 

1 ^ r fp(ri)]o„ [p(r3)]„i [p(r2)]w 



(4) 
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I [p(r2)]o„ [p(ri)]„Jp(r3)];o ^ [p(r3)]o^ [p(r2)]„; [p(ri)]^o 
[t^' - Uno - irj] [u" + - \rj\ [u" - u^o - ii]] [u + ujiq + ir]] 

I [P(ri)]0n [P(r2)]„; [p(r3)];o ^ [p i^-j)] On [P l)] nl iPi^ 2)] ip 

[u - Uno + ir]] [^" + ^10 - iv] [^" - ^no - iv] + ^lo - 'M 

^ W2)]on [Piy?)]nl Wl)]lO 1 _ /gx 

[u' - ujno - iri] [u + + ir/] J ■ 

Here, Vt is the normalization volume, p(r) is the particle-density operator, and 
[p(r)]^; are matrix elements taken between the exact many-electron states \E'„ 
and '^i of energy En and Ei. $o and Eq represent the exact many-electron 
ground state and energy, respectively, Uni = En — Ei, u" = —{u + u'), and rj 
is a positive infinitesimal. 

The total energy lost by the probe particle is simply 

which can also be obtained by first considering the probability Pq for the 
probe particle to transfer momentum q and energy uj to the many-electron 
system and then multiplying Pq^ ^^ by the energy transfer u and summing over 

all q and a; 123113 nails!. 

The many-body ground and excited states of a many-electron system 
are unknown; hence, the exact linear and quadratic density-response func- 
tions are difficult to calculate. In the framework of time-dependent density 
functional theory [TDDFT] |in], the exact density-response functions are 
obtained from the knowledge of their noninteracting counterparts and the 
exchange-correlation (xc) kernel /xc(r, r', u) which equals the second func- 
tional derivative of the unknown xc energy functional i^^cl?^]- In the so-called 
time-dependent Hartree approximation or RPA, the xc kernel is simply taken 
to be zero. 

In the case of a noninteracting Fermi gas, the ground state is obtained 
by simply filling all the single-particle states of noninteracting electrons be- 
low the Fermi level. When acting on the ground state, the particle-density 
operator p(r) produces single-particle transitions in which a given particle is 
scattered from some state inside the Fermi sea to a state outside. Hence, the 
linear and quadratic density-response functions of Eqs. (jH) and Q reduce to 
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their noninteracting counterparts: 



0.(ri)0*(ri)0,(r2)0:(r2) 0.(r2)0*(r2)0,(ri)0*(ri^ 



uj — Uji + ir] 



UJ + UJji + ITj 



(7) 



and 



r°(ri, rs, rg; u;') = ^ /, 



0,(ri)0*(ri)0,(r3)0I,(r3)0fc(r2)0:(r2) 
[u; - Uji + ir]] [cu' + ujki - \rf\ 



^ 0,(r2)0-(r2)0,(ri)0^(ri)0fe(r3)0:(r3) 0,(r3)0*(r3)0,(r2)0Ur2)0fc(ri)0*(ri; 
H r-T"^ . , , .. . — \ 



[uj' - ujji - irj] [u" + uJki - irj] 



" - iOji - irj] [u + UJki + i^] 



^ 0.(ri)0-:(rO0,(r2)0^(r2)0,(r3)0:(r3) 0.(r3)0*(r3)0,(ri)0^(ri)0fc(r2)0:(r2) 
H — r-Ti—r. h - 



+ 



[uj - Uji + irf\ [uj" + UJki - yq] 
0,(r2)0;(r2)0,(r3)0^(r3)0fc(ri)0*(rO 
[uj' - ujji - irj] [uj + UJki + i^?] 



[u;" - c^jj - ir/] [c^' + cjji - irf\ 



(8) 



where cUjj = Ej — and fi are Fermi-Dirac occupation factors. At zero tem- 
perature fi = Q{eF — £i), Sf being the Fermi energy, and 6(x), the Heav- 
iside step function. The single-particle states 0j(r) and energies Si entering 
Eqs. ([7j) and (jH)) are usually chosen to be the eigenfunctions and eigenval- 
ues of an effective Hartree |17j, Kohn-Sham |3H|, or quasiparticle |49i Eni 
hamiltonian 

In the RPA, the induced potential V^^'^iv, t) is simply obtained as the 
potential induced in a noninteracting Fermi system by both the probe-particle 
charge density p^^'*(r, t) and the induced electron density p^^'^{r,t), i.e., 



I dr' I dt' I "^e-^-(*-*') / rfri / dv^v{v,v^ 

J J-oo J-oo ZTT J J 

°° duj' 



2tt 



X |x°(ri,r3;a;) + J dr" J^^ dt" J_ 
X /.r.y0(r,r„r3;.,.0.(r„0 

Assuming that 

V^^"'^(r,t) = J dr'v{r,r')p'''\r',t), 



i{u+ui'){t'-t") 



(9) 
(10) 
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one easily finds that the induced potential V^^'^{r, t) of Eq. IQ is of the form 
of Eq. (j2)) with the exact interacting linear and quadratic density-response 
functions being replaced by the following integral equations: 

X(r. .) = x»(., r'; .) ^ j i., j i.-.x>, (11) 
and 

y(r. = /dr,/<ir,/<ir3/f(r,.,.)yO(r.,r..r3;.,.') 

X K(r2,r';-u')K(r3,i";u: + ui'), (12) 
where K{r, r'; a;) is the so-called inverse dielectric function 

K{r, r'; to) = 5{r - r') + J dvi x(r, ri; cu) v{ri, r'). (13) 
2.1 Uniform electron gas 

In the case of a uniform FEG, there is translational invariance in all direc- 
tions. Hence, Eq. ^ yields 

where L is the normalization length, and {—dE/dx) is the energy loss per 
unit path length of the projectile, i.e., the so-called stopping power of the 
electron system: 

dE ^ f dq du ^ , , 

— ImKq^^ -\- 2n Zi / — - — — / Iml^ (^._q^ V^^ ^q-qi 

5 (cu? - qi 



Vq = 47r/g^, Kf^,J, and Fq^,tji;q2,w2 being Fourier transforms of the bare 
Coulomb interaction v{r,r'), the inverse dielectric function K{r,r',uj), and 
the quadratic density-response function F(r, r', r"; tu, a;'), respectively. 
In the RPA, 

^q,a; = 1 + Xq.o; f^q, (16) 
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and 

^l,q2;t<Jl,a;2 ~ -^qi,a;i ^i,q2;LJi,a;2 -^-qi,-i^l -^qi +q2 ,1^1 +0^2 5 (l^) 

Xq^^ and ^qj^qjjoji.i^a being Fourier transforms of the noninteracting linear and 
quadratic density- response functions x^lri, 1*2; and F°(ri, r2, ra; cj, tu') [see 
Eqs. (jlj) and (jH))]- Noting that for a uniform electron gas the single-particle 
states 0i(r) entering Eqs. ((Zj) and (jHl) are simply plane waves of the form 

0k(r) = -^e^*^- (19) 



analytic expressions for both x^ ,^ and ^qi,q2;a;i,(^2 '^^'^ obtained. Xq,a; 
the Lindhard function jHlllZj. Explicit expressions for the real and imaginary 
parts of the noninteracting quadratic density-response function Y^^ ,^^.^_^^^^ 
where reported in Refs. and [121 123, respectively, and an extension to 
imaginary frequencies was later reported in Ref. 

2.1.1 High- velocity limit 

For high projectile velocities, the zero-point motion of the electron gas can 
be neglected and it can be considered, therefore, as if it were at rest. Thus, 
in this approximation all energies Si entering Eqs. ((Tj) and ^ can be set 
equal to zero. If one further assumes that v"^ » Up, uop being the plasmon 
frequency for which i^o.i^ diverges, Eq. (fT3j) is found to yield [22] 

dE 4:7mZ^ , , , , 

- T. " + P"' 

where 

U = \^— (21) 

Up 

and 



Li = 1.42^ In , 22) 



n being the density of free electrons in the target. The first term, which 
for a relatively low projectile charge Zi gives the main contribution to the 
stopping power, is proportional to Zl and has the same form as the Bethe 
formula for the inelastic stopping power of atoms |^ as long as the plasma 
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frequency Up is replaced by the mean excitation energy of electrons in the 
atom. The second term, which originates the so-called Barkas effect, i.e., the 
dependence of the stopping power on the sign of the projectile charge, has 
been found to yield excellent agreement with stopping-power measurements 
at high velocities [T^ 



2.1.2 Double-plasmon excitation 

In the RPA, the linear- response contribution to the stopping power of Eq. (jl5p . 
which is proportional to Z^, is fully originated in the creation of single e-h 
pairs and plasmons. Furthermore, the contribution to the actual (beyond 
RPA) Zl stopping power that is due to coherent multiple excitations such 
as double plasmons is expected to be relatively small. Nevertheless, accu- 
rate measurements of electron energy-loss spectra showed evidence for the 
existence of coherent double-plasmon excitations [T71 ITH]. 

In the framework of many-body perturbation theory, one first defines the 
scattering matrix as a time-ordered exponential in terms of the perturbing 
hamiltonian and field operators |47|. Then, one considers the matrix elements 
corresponding to the process in which the recoiless probe particle carries 
the system either from an initial state a^$o to a final state a/$o (single 
excitation) or from an initial state a^a^$o to a final state Oj^aj^^o (double 
excitation): 



and 





\af 


Saf\ 


|$o > 




\S\% > 



su = (23) 



< %\af,af,Sala+\% > 
- < $o|5|$o > ' ^^^^ 

where Oj and af are annihilation and creation operators, respectively, and $o 
represents the vacuum state. Finally, one calculates the probability for the 
probe particle to transfer momentum q and energy uj to the many-electron 
system by moving either one, two, or more particles from inside the Fermi 
sea to outside: 



p^,^ = 2E/sE(i-/p)l%l'C-. + 4EEE/siE/s.E(i-/pJ 

s p qi a;i si S2 pi 

^ E(l ~ /pa) \Sflj2;iui2f ^ti,Pi-si^t-quP2-s2 ' (25) 



P2 
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^q,q' ~ ^q,q''^'^,'^' being the symmetric Kronecker S symbol. If the probe- 
particle is not a heavy particle, energy conservation should be ensured by 
introducing recoil into the argument of the delta functions. If the probe is an 
electron, a step function should also be introduced to ensure that the probe 
electron does not lose enough energy to fall below the Fermi level. 

The stopping power of the electron system is obtained by first multiplying 
the probability Pq^^ of Eq. (j25j) by the energy transfer lj and then summing 
over all q and u: 

-f "lEE-^",."- (26) 

A careful analysis of the various contributions to the probability Pq^^^ that 
are proportional to Zf and Zf yields, to third order in the RPA screened 
interaction v^K^^,^, an expression for the stopping power that exactly coin- 
cides with Eq. (jl5|) [44]. Alternatively, the inverse mean free path is obtained 
as follows 

A-^ = 7EE^q,.. (27) 

The lowest-order Zf contribution to the probability Pq^ for a probe- 
electron {Zi = —1) to excite a double plasmon, which is of fourth order in 
the RPA screened interaction v^K^^^, is found to be given by the following 
expression: 



= ^^^Z'^v^^K^lY.Vn-m I ^ ImfsTq,,^, ImAq_q,,^_^, 



X |yq,q,;^,^J^5(g°-/ + CUv-q)e(c^v-q-eF), (28) 

where = k'^/2. If one approximates both the linear and quadratic density- 
response functions entering Eq. by their low-q limits, one only keeps the 
high-velocity limit of this probability in an expansion in terms of the inverse 
velocity, and one introduces this limit into Eq. ()27|) . one obtains the following 
result for the Zf contribution to the inverse mean free path that is due to 
the excitation of a double plasmon PUIIH^ : 



Vs being the so-called electron-density parameter rs defined by the relation 
1/n = (47r/3)(r5 Oq)^. Numerical study shows [21] that introduction of the 
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full RPA linear and quadratic density-response functions into Eq. (j28|) yields 
a result for X^p which has, in the high- velocity limit, the same velocity de- 
pendence as the approximation of Eq. (P^ . though for = 2.07 the full 
RPA X^} is found to be larger than this approximation by a factor of 2.16. 



2.2 Bounded electron gas 

In the case of charged particles moving inside a solid, nonlinear effects are 
known to be crucial in the interpretation of energy-loss measurements. How- 
ever, these corrections have been shown to be less important when the 
charged particle moves outside the solid IHE!- Hence, in the case of a bounded 
three-dimensional electron gas we restrict the calculations to linear-response 
theory. Assuming translational invariance in two directions, which we take 
to be normal to the z axis, to first order in the external perturbation (linear- 
response theory) the energy loss of Eq. (0) may be expressed as follows [SZl 

-AE = -^/-^/ dt dt' duu 



n J (27r)2y-oo 

X e-'("-'it|-"ii)(*-*')lmiy[z(t),2(t');q||,cu], (30) 

where qy and vy are components of the momentum transfer and the velocity 
in the plane of the surface, z(t) represents the position of the projectile 
relative to the surface, and 

W{z,z';qii,uj) = v{z, z\ci\\)+j dzi j dz2v{z,Zi\(i\\)x{zi,Z2](\\\,uj)v{z2,z'\q,\i), 

(31) 

v{z,z']C{\\) = (27r/g||)e ''n ^ ' and x(2;, 2;'; qy, cu) being two-dimensional Fourier 
transforms of the bare Coulomb interaction and the density-response func- 
tion, respectively. 
In the RPA, 

X(...';q,..) = X°(...';q„.-) + /<i../d.«V...;q|,.-) 

X t;(zi,Z2;q||)x(2;2,2;';q||,u;), (32) 

where 

/(.,.';q,„.) ^ 2 E0.W0K^)0.<^O0*(/) / (^^^^TFT^- 

(33) 
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Here, 

Ei = ei + ^ (34) 

and 

E, = + &±^, (35) 

the wave functions (piiz) and energies Si, which describe motion normal to 
the surface, being the eigenfunctions and eigenvalues of a one- dimensional 
Hartree, Kohn-Sham, or quasiparticle hamiltonian. 

Eq. (jHUj) gives the energy that a charged particle moving with constant 
velocity along an arbitrary trajectory loses due to electronic excitations in an 
electron system that is translationally invariant in two directions, as occurs 
in the case of a simple metal surface modeled by jellium. 



2.2.1 Parallel trajectory 

In the glancing incidence geometry ions penetrate into the material, they 
skim the outermost layer of the solid, and are then repelled by a repulsive, 
screened Coulomb potential, as discussed by Gemmell jSH]- Through use of 
the appropriate effective potentials the ion trajectory z{t) can be calculated 
and the energy loss is then obtained from Eq. (|3U|1 . Under extreme grazing- 
incidence conditions, incident charged particles can be assumed to move with 
constant velocity v along a definite trajectory at a fixed distance z from a 
jellium surface. Eq. (jHUI) then yields Eq. (Hlj), as in the case of charged 
particles moving in a uniform electron gas, but with the energy loss per unit 
path length now being given by the following expression |57j : 

~ Ih ^ ~v^' J J2^ Jo duJUjlmW{z,z;q\\,uj)5{uj-qwyr). (36) 

2.2.2 High- velocity limit 

At high velocities, the energy-loss spectrum for charged particles moving out- 
side a solid is known to be dominated by long-wavelength (gy — > 0) surface- 
plasmon excitations [HH]- In this long- wavelength limit, the imaginary part 
of Eq. dSH) yields jBD] 

lmW{z,z';qn,uj) = e"^" a;, 5(a; - u;,), (37) 

^11 
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where Ug = ujp/\/2, ujp = {Ann)^^^ being the classical plasma frequency of a 
uniform electron gas of density n. 

Introducing Eq. (jHTj) into Eq. (jHBj) . one easily reproduces the classical 
expression of Echenique and Pendry : 



dE „o u 



2 



dx Zt:;^^oi2u;sz/v), (38) 

where Kq is the zero-order modified Bessel function. For large values of z 
{z » v/ujs), Eq. reduces to 

-f-^? I; (39) 
2.3 Periodic crystals 

For a peridic crystal, we introduce the following Fourier expansion of the 
linear density-response function: 

1 BZ 

q G,G' 

where the first sum runs over q vectors within the first Brillouin zone (BZ), 
and G and G' are reciprocal lattice vectors. 
In the RPA, 

XG,G'(q>^) = XG,G'(q>^) + EXG,G"(q>^)^q+G"XG",G'(q,w), (41) 

G" 

where fq+G and Xg ci^^^) represent the Fourier coefficients of the bare 
Coulomb interaction v{r, r') and the noninteracting density-response function 
X^{r, r'; u), respectively: 

4:71 



and 



BZ 

XG,G'(q5 ^) = ^ E E (/k,n - /k+q,nO 
k n,n' 



Ek,n - £^k+q,n' + h{uj + IT]) 
14 



0k,n and i?k,n being Bloch eigenfunctions and eigenvalues of a three-dimensional 
Hartree, Kohn-Sliam, or quasiparticle hamiltonian. 

The stopping power of a periodic crystal is obtained by first introducing 
Eqs. (0) and © into Eq. ((T)), and then introducing Eq. (jT)) into Eq. 
Within linear-response theory, i.e., to first order in the external perturbation, 
the result is of the form of Eq. (jl4p with the energy loss per unit path length 
being given by the following expression: 



J pi o 72 BZ 

- ^ EEE^e'^'^-'^.+G+G' ImirG,G+G'Jq,cu), (44) 



dx ffi qGG 



where b is the impact vector of the projectile, the sum J2g'j_ is restricted 
to those reciprocal-lattice vectors that are perpendicular to the projectile 
velocity (G'j_ ■ v = 0), = (q+ G) ■ v, and KG G>{q, u) represent the Fourier 
coefficients of the inverse dielectric function K{r, r'; uj): 

Kg,g' (q, t^) = 5g,g' + Wq+G Xg,g' (q, c^) • (45) 

The stopping power of Eq. (jl^ is the so-called position-dependent stop- 
ping power. The most important contribution to this quantity is provided 
by the term = 0, the magnitude of the other terms depending on the 
direction of the velocity. For a few highly symmetric or channeling direc- 
tions non-negligible corrections to the G'j^ = contribution are found, thus 
exhibiting the characteristic anisotropy of the position-dependent stopping 
power. However, for those directions for which the condition of G'j^ ■ v = 
is never satisfied only the G^ = term contributes and one finds 

2 Z"^ 

= ^ E E ^ ^q+G Im i^cG (q, w) • (46) 

random Q G 

This is the so-called random stopping power, which is also obtained as the 
average over impact parameters of the position-dependent stopping power 
of Eq. For simple metals like Al the diagonal elements of the inverse 

dielectric matrix Kci^G{q^,uj) are rather isotropic, in which case there is little 
dependence of the random stopping power on the direction of the projectile 
velocity. 



dE 
dx 
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3 Results 



In this section, we review existing calculations of the stopping power of both 
an infinite (uniform) and a semi-infinite (bounded) free gas of interacting 
electrons, double-plasmon inverse mean free paths, and the linear (lowest- 
order) first-principles stopping power of real Al and Si. New calculations 
of the various contributions to the stopping power due to the excitation of 
e-h pairs and plasmons are also reported, existing double-plasmon calcula- 
tions are extended to low velocities where both recoil and the probe-particle 
statistics play a role, and existing first-principles calculations of the stopping 
power of Al and Si are extended by allowing transitions with large values of 
the momentum transfer. 

3.1 Uniform electron gas 

Here we consider a uniform gas of interacting electrons, the electron density n 
being equal to the average electron density of valence electrons in aluminum 
metal (r^ = 2.07), for which the Fermi momentum [qp = (Stt^tt,)^/'^] and 
bulk plasma frequency [tOp = (4?™)^/^] are qp = 0.927 and tOp = 15.8 
eV, respectively. We set Zi = 1 (probe bare protons) for the stopping- 
power calculations and Zi = —1 (probe electrons) for the calculations of 
the double-plasmon inverse mean free path. Our results can then be easily 
extended to arbitrary values of Zi, as linear and quadratic contributions to 
the stopping power and the inverse mean free path are proportional to Zf 
and Zf, respectively. 

3.1.1 Stopping power 

The second-order stopping power of a homogeneous FEG is given by Eq. (|15p. 
In Fig. [T] we show, as a function of the projectile velocity, our full RPA 
calculations for the separate Zf and Zf contributions to the stopping power 
of Eq. (fT3j) . We note that both Zf and Zf contributions to the stopping 
power exhibit a linear dependence on the projectile velocity up to velocities 
approaching the stopping maximum. This linear dependence is also exhibited 
by full nonlinear DFT calculations of the stopping power of a FEG j2S] 
and by recent measurements of the electronic energy loss of protons and 
antiprotons [T^ . 

The linear dependence of the Zf contribution to the stopping power is a 
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V (a.u.) 



Figure 1: Solid lines represent RPA Zl and contributions to the stopping 
power of Eq. (fl3|). for Zi = 1 and = 2.07, as a function of the velocity of 
the projectile. Dashed and dashed-dotted lines represent the RPA [—dE/dx]^ 
and [—dE/dx]^ contributions to the Zf stopping power, respectively. The 
dotted line is the high- velocity limit dictated by the second term of Eq. (|^. 
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consequence of two competing effects. On the one hand, there is the effect 
of one-step single excitations, hke those entering the Z\ contribution to the 
stopping power, but now generated by the quadratically screened potential 
of the probe particle. On the other hand, there is the effect of second-order 
two-step single excitations generated by the linearly screened probe potential. 
In the RPA, these contributions to the Z\ effect are given by the following 
expressions PU] : 



dE 
dx 



X Re Vq™_q^ V^^ Re /STq,,^, t;q_q, Re irq_qi,^_c.i5 [^l - " V 



(47) 



and 



dE_ 
dx 



V J (zvr)-* JO In J (/tt)* j-oo ztt 



X -^q,a;; — qi,— ciJi ^qi Re-K'qxjtJi ^q— qi Re -^^q— qj^ 



6 {uj'( - qi ■ v) , 



(48) 



where ^(^5!qi,wi is the time-ordered counterpart of the retarded non-interacting 
quadratic density-response function Y^^^.^_^^^_^, and 



H. 



q,a;;qi,i^i 



27rsgn(co')P 



d^k 



(2 



/k 



LlJi + UJ]i — a;k+qi 



+ 



^^k + t^k+q2 

(49 

with q2 = q — qi and 002 = w — Ui. Double excitations also contribute, within 
RPA, to the Zf stopping power of Eq. (fT^. but they are found to be zero in 
the low and high velocity limits and small at intermediate velocities. 

The RPA [—dE/dx]^ and [—dE/dx]'^ contributions to the stopping power 
of a PEG are represented in Fig.[T]by dashed and dashed-dotted lines, respec- 
tively. The [—dE/dx]'^ contribution from losses to two-step single excitations 
(dashed-dotte line) is very small at high projectile velocities where the veloc- 
ity distribution of target electrons can be neglected. At high velocities, the 
surviving [—dE/dx]^ contribution approaches the total Zf stopping power, 
which is very well reproduced by the second term of Eq. (j^UI) (dotted line of 
Fig.lH). 
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V (a.u.) 



Figure 2: The solid line represents the RPA Zl (linear) contribution to 
the stopping power of Eq. (fT3j). for Zi = 1 and = 2.07, as a function 
of the velocity of the probe particle. Dashed and dotted lines represent 
contributions from the excitation of plasmons and e-h pairs, respectively. 



Now we focus on the role that the excitation of e-h pairs and plasmons 
plays in the energy loss process. First, we restrict our calculations to the 
Zl (linear) contribution to the stopping power of a FEG. Fig. |21 exhibits the 
separate RPA e-h pair and plasmon contributions to the Z\ stopping power 
of Eq. (jl5j) . This figure shows that the contribution from losses to plasmons 
is smaller for all projectile velocities than the contribution from losses to e-h 
pairs, which is especially true at high electron densities, although both con- 
tributions coincide in the high-velocity limit. This is the equipartition rule, 
which appears straightforwardly in the static-electron gas approximation. In 
this approximation, plasmon and e-h pair contributions to the energy loss 
are typically separated according to whether the momentum transfer is be- 
low (g < Qc) or above [q > qc) the critical momentum qc where the plasmon 
dispersion enters the e-h pair excitation spectrum. For an electron gas not 
at rest, the equipartition rule was formulated by Lindhard and Winther 

For a momentum transfer that is smaller than q^, both plasmon and e-h 
pair excitations contribute to the full RPA energy loss, though contributions 
from losses due to the excitation of e-h pairs are very small. For q > q^ how- 
ever, only e-h pair excitations contribute. This is illustrated in Fig. where 
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Figure 3: The solid line represents the RPA Zl (linear) contribution to 
the stopping power of Eq. (fT3jl . for Zi = 1 and = 2.07, as a function 
of the velocity of the probe particle. Dashed and dotted lines represent 
contributions from momentum transfers below (g < Qc) and above (g > Qc) 
the critical momentum Qc where the plasmon dispersion enters the e-h pair 
continuum. 
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the total Zl stopping power is separated according to whether losses corre- 
spond to momentum transfers below (g < q^) or above (g > q^) the critical 
momentum q^. In Fig. |3J the total Zl stopping power is separated accord- 
ing to whether losses come from momentum transfers below (g < ^/2ujp) or 
above (g > a/2ca^) the critical momentum ^/2uj^, which is the low-density 
limit of qc- In this case, there is exact equipartition for all velocities above 
the stopping maximum. This equipartition is also found to be exact in the 
high- velocity limit, by using Coulomb scattering of independent electrons 
with gmin = ojpl V or by assuming that independent electrons are scatterd by 
a velocity-dependent Yukawa potential with screening length proportional to 

Up/V. 

While the RPA [—dE/dx]'^ contribution to the Zf stopping power [see 
Eq. fl48|) ] is entirely due to the excitation of e-h pairs, the [—dE/dx]^ contri- 
bution to the zf stopping power is originated in the excitation of both e-h 
pairs and plasmons. Hence, we have split [—dE/dx]^ into these contributions 
and have found the result shown in Fig. El by dotted and short-dashed lines, 
respectively. This figure shows that contributions to the Zf effect coming 
from losses to plasmons is relatively smaller than in the case of the Zf term 
(see Fig. 12)), especially at high velocities. Hence, collective excitations ap- 
pear to be reasonably well described with the use of linearly screened projec- 
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Figure 5: The solid line represents the RPA Zf (quadratic) contribution to 
the stopping power of Eq. (fT3jl . for Z\ = \ and = 2.07, as a function of 
the velocity of the probe particle. Long-dashed and dashed-dotted lines rep- 
resent the RPA \—dEjdx^ and \—dEjdx^ contributions to the Z\ stopping 
power, respectively. Short-dashed and dotted lines represent contributions 
to \—dE I dx^ from the excitation of plasmons and e-h pairs, respectively. 



tile potentials. The equipartition rule, valid within first-order perturbation 
(linear-response) theory cannot be extended to higher orders in the external 
perturbation. 

In order to account approximately for the Z\ effect coming from both 
the conduction band and the inner shells, a local-plasma approximation was 
used in Ref . [12] , by assuming that a local Fermi energy can be attributed to 
each element of the solid. The experimental differences between the stopping 
power of silicon for high-velocity protons and antiprotons ^2] were succes- 
fuUy explained in this way. Proton and antiproton stopping powers of a 
variety of solid targets have been reported recently for projectile velocities 
near and below the stopping maximum A comparison of our theory 

with these experiments requires the inclusion of charge-exchange processes, 
xc effects, and higher-order nonlinear terms. Work in this direction is now 
in progress 
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Figure 6: RPA double-plasmon inverse mean free paths of electrons (solid 
line) and positrons (dashed line) for = 2.07, versus the velocity of the 
projectile, as obtained from Eqs. and by either including (electrons) 
or excluding (positrons) the step function 6(c<Jv-q — Sp)- The dotted line is 
the high- velocity limit dictated by Eq. (jSUj) . 



3.1.2 Double-plasmon excitation 

Single-plasmon contributions to the electron inelastic mean free path of swift 
electrons have been calculated for many years, both in the high-velocity 
limit jni] and in the full RPA jHS]. Calculations of the double-plasmon 
contribution to the electron inverse mean free path have been reported in 
Refs. |19j.20ti25. Fig.inishows our full RPA calculation of the double-plasmon 
inverse mean free path of swift electrons (solid line) interacting with a FEG, 
as obtained after introducing Eq. ()28|) into Eq. (jUj). Also shown in this fig- 
ure is the double-plasmon inverse mean free path of swift positrons (dashed 
line), as obtained by simply removing the step function of Eq. fl28p . and the 
high- velocity limit dictated by Eq. (^^1) multiplied by a factor of 2.16 (dotted 
line), i.e.. 



At high velocities of the projectile the zero-point motion of the target can 
be neglected. Hence, at these velocities the effect of the Pauli restriction, 
which only applies to electron probes, is negligible, and the behaviour of 




(50) 
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the double-plasmon inverse mean free path is independent of the particle 
statistics. On the other hand, it is interesting to notice that the high- velocity 
formula of Eq. (j^Uj) gives an excellent account of the full RPA result for both 
probe electrons and positrons in a wide range of projectile velocities. In 
particular, for Al and a probe-electron energy of 40keV Eq. (j^UI) yields a 
ratio for the double relative to the single plasmon inverse mean free path of 
1.9 X 10~^, in agreement with experiment |18j . 

3.2 Bounded electron gas 

First, we consider a jellium slab of thickness a normal to the z axis, consisting 
of a fixed uniform positive background of density 

{n, —a < z <{) 
(51) 
0, elsewhere, 

plus a neutralizing cloud of interacting electrons of density n{z). The positive- 
background charge density n is expressed in terms of the Wigner radius 
[l/n = (47r/3)(rs ao)^], which we take to be = 2.07. 

To compute the interacting RPA density- response function of Eq. (jH^ . we 
follow the method described in Ref. [nE]- We first assume that n{z) vanishes 
at a distance Zq from either jellium edge [HZl, and expand the wave functions 
(j)i{z) in a Fourier sine series. We then introduce a double-cosine Fourier 
representation for the density-response function, and find explicit expressions 
for the stopping power of Eq. ()36|) in terms of the Fourier coefficients of 
the density-response function jHZ|. We take the wave functions (pii^z) to be 
the eigenfunctions of a one-dimensional LDA hamiltonian with use of the 
Perdew-Zunger parametrization jnS] of the Quantum Monte Carlo xc energy 
of a uniform FEG j69i 

Finally, the stopping power of a semi-infinite FEG is obtained with the 
use of the following relation: 

dE [-dE/dx] (a^) + \-dE/dx\ (a„) + [-dE/dx] (a+) 
" = 3 ^^^^ 

where a„ is the threshold width for which the nth subband for the z motion 
is first occupied, a~ = a„ — Air/4, and a+ = a„ + Ai?/4, Xp = 27r / {Sn'^nY^^ 
being the Fermi wavelength [20] • Converged results have been found with 
the use of slabs with n = 12, for which a ~ 5 — 6Ai?. 
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z (a.u.) 

Figure 7: The solid line represents the RPA stopping power of Eq. fj36p for 
Zx = 1, = 2.07, and f = 2fo, as a function of the distance z between the 
surface and the projectile. The solid is in the region z < 0. The dashed line 
represents the LDA stopping power, as obtained by assuming that the actual 
stopping power of Eq. (jHUI) can be approximated by that of a uniform PEG 
with the local density ni^z). 
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Figure 8: The solid line represents the RPA stopping power of Eq. flHfjj) for 
Zi = 1, = 2.07, and z = 3ao ~ ^f/'2, as a function of the projectile 
velocity. The dashed line represents the LDA stopping power, as obtained 
by assuming that the actual stopping power of Eq. can be approximated 
by that of a uniform PEG with the local density n{z). The dotted line is the 
high- velocity limit dictated by Eq. (j38j) . 



Fig. m depicts our full RPA calculation of the stopping power of a semi- 
infinite PEG, as obtained from Eq. (jHUjl for protons {Zi = 1) moving with 
speed V = 2vq parallel to the surface. In the interior of the solid, where the 
electron density is taken to be constant, the stopping power coincides with 
the Zf (linear) RPA stopping power of a uniform PEG (see Fig. H}. Outside 
the solid, the stopping power decreases with the distance z between the 
surface and the probe-particle trajectory. Also plotted in Pig. [7|is the result 
of assuming that the stopping power for a charged particle that moves at a 
distance z from the surface can be approximated by that of a uniform electron 
gas with the local density n{z) Pig. [7| clearly shows that this often-used 
local-density approximation yields an inaccurate description of the position- 
dependent stopping power, due to the intrinsic nature of surface-induced 
single and collective excitations not present within this approach. 

As the velocity increases, the energy-loss spectrum of charged particles 
moving far from the surface into the vacuum is dominated by long-wavelength 
excitations. In this limit, Eq. yields the classical stopping power dictated 
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by Eq. (jHS|) . This is illustrated in Fig. |S1 where the velocity-dependent RPA 
stopping power of Eq. is represented together with the classical result [see 
Eq. (jHHj) ] for a proton {Zi = 1) moving in the vacuum at a constant distance 
z = 3 ao ^ Xp/2 from the surface. At low velocities, the energy- loss spectrum 
is dominated by intermediate and short wavelength excitations, even far from 
the surface into the vacuum, and a combination of the actual electronic sel- 
vage at the surface with the intrinsic nature of surface-induced excitations 
increases the energy loss with respect to that predicted by Eq. (|HHj). At high 
velocities, the energy-loss spectrum is dominated by surface-plasmon exci- 
tations and the full RPA stopping power nicely converges with the classical 
result. As in Fig. [TJ the local-density approximation is also represented in 
this figure, showing that this often-used approximation cannot account for 
the energy loss originated in surface-induced excitations, not even at low 
velocities where the energy loss is entirely due to the excitation of e-h pairs. 

3.3 Periodic crystals 

In the case of periodic crystals, we first expand the one-electron 0k,n(r) eigen- 
functions in a plane-wave basis, 

0k,n(r) = ^E^M(G)e'(''+«)-, (53) 

with a kinetic-energy cutoff that varies from 12 Ry in the case of Al (~ 100 G- 
vectors) to 16 Ry in the case of Si (~ 300 G- vectors). The coefficients Mk,n are 
evaluated by solving the Kohn-Sham equation of DFT in the LDA with use 
of the Perdew-Zunger parametrization [68, of the Quantum Monte Carlo xc 
energy of a uniform FEG [Hni • The electron- ion interaction is based on the use 
of an ab initio non-local, norm- conserving ionic pseudopotential [^. Then, 
Eq. (j43|) is used to evaluate the Fourier coefficients of the noninteracting 
density-response function, and a matrix equation [Eq. (|4ip] is solved for the 
Fourier coefficients of the interacting RPA density-response function. 

Here we present new calculations for the random stopping power of Al and 
Si, which represent an extension of existing first-principles calculations [3H1 
HO] by allowing transitions with larger values of the momentum transfer. 
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Figure 9: First-principles RPA calculation of the random stopping power 
of valence electrons in Al (solid circles) and Si (open circles) for protons 
and antiprotons {Zf = 1), versus the projectile velocity, as obtained from 
Eq. ()4fi|l . These results have been found to be rather insensitive to the choice 
of the direction of the projectile velocity. The solid line represents the Zf 
(linear) stopping power of a uniform FEG with = 2. 



3.3.1 Random stopping power 

Fig. IHl shows our first-principles first-order (linear-response) RPA calculation 
of the random stopping power of valence electrons in Al (solid circles) and 
Si (open circles) for protons and antiprotons {Zf = 1), as obtained from 
Eq. [75] . Since the electron density of valence electrons in Al (r^ = 2.07) 
and Si (r^ = 2.01) is nearly the same, within a FEG model of the solid the 
stopping powers of Al and Si are expected to coincide (solid line of Fig. IH}. 
Nevertheless, our first-principles calculations indicate that this is not the 
case. At low projectile velocities (where only valence electrons contribute 
to the energy loss of the projectile), the stopping power of Si is consider- 
ably smaller than that of Al, in agreement with available measurements of 
the stopping power of this materials for either protons |73] and antipro- 
tons [13j. While the stopping power of Al is found to be slightly larger than 
in the case of a FEG, the band gap of Si yields a stopping power of this 
material that is smaller than in the case of a FEG. 

At high velocities, well above the stopping maximum, the sum over the 
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frequency uj in Eq. (j46|) can be replaced by an integration over all positive 
frequencies, and the sum rule 



/■oo 

doj uj\m.K^Q^i{c{,uj) = 27r^?2G-G' (54) 



Jo 



(nc represents the Fourier components of the density, ^0=0 being the average 
electron density n of the crystal) yields a stopping power which depends on 
n but not on the details of the band-structure of the target material: 



dE 
dx 



— ^nln— . (55) 



- random ^ 



Hence, at high velocities the stopping powers of valence electrons in Al and 
Si both coincide with that of a FEG with the same electron density (see 
Fig.inD. 

While at low velocities the contribution to the total energy loss due to 
excitation of inner-shell electrons is negligible, at velocities larger than the 
Fermi velocity it is necessary to allow for this contribution. The cross sections 
for the ionization of inner shells in Al were obtained by Ashley et al. [70] in 
the first-Born approximation utilizing atomic generalized oscillator strength 
functions. By adding the contribution from core electrons to that of valence 
electrons (this contribution was calculated within a FEG model of the solid) 
these authors found a nice agreement with experiment. Good agreement 
with experiment was also shown in Ref. |2j by adding to the valence-electron 
energy loss of a FEG with = 2.01 the energy loss from core electrons in Si 
as taken from Walske's calculations [77j . 



4 Summary and conclusions 

We have presented a survey of current investigations of various aspects of 
the interaction of charged particles with solids. 

In the framework of many-body perturbation theory, we have studied 
the nonlinear interaction of charged particles with a free gas of interacting 
electrons. We have presented general procedures to calculate the nonlinear 
potential induced by charged particles moving in an inhomogeneous electron 
system, the Zf contribution to the stopping power of a FEG, and double- 
plasmon excitation probabilities. 
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Our calculations for the RPA stopping power of a FEG indicate that 
for velocities smaller than the Fermi velocity the stopping power is, up to 
third order in the projectile charge, a linear function of the velocity of the 
projectile. Our calculations also indicate that the high-velocity limit dictated 
by Eq. ()2U|) gives an excellent account of the full RPA result in a wide range of 
projectile velocities. By assuming that a local Fermi energy can be attributed 
to each element of the solid target, the experimental differences between the 
stopping power of silicon for high- velocity protons and antiprotons fT2] were 
succesfuUy explained in Ref. P^ . 

New calculations of the various contributions to the second-order stop- 
ping power of a uniform FEG coming from the excitation of e-h pairs and 
plasmons have been reported. We have found that the equipartition rule, 
valid within first-order perturbation (linear-response) theory, cannot be ex- 
tended to higher orders in the external perturbation. We have also found 
that contributions from collective excitations to the Zf term are small. 

New RPA double-plasmon inverse mean free paths of electrons and positrons 
have been evaluated, with explicit introduction of recoil and probe-particle 
statistics. The high-velocity limit of Eq. (j5Up is found to give an excellent 
account of the full RPA double-plasmon inverse mean free paths in a wide 
range of projectile velocities. This formula [Eq. ()50|)] yields for Al and a 
probe-electron energy of 40 keV a ratio for the double relative to the single 
plasmon inverse mean free path of 1.9 x 10~^, in agreement with experi- 
ment [TH] . 

We have reviewed existing self-consistent calculations of the energy loss 
of charged recoiless particles moving parallel to a plane-bounded FEG, in 
the framework of linear response theory. In the high-velocity limit and for 
charged particles moving far from the surface into the vacuum the actual stop- 
ping power is found to converge with the classical limit dictated by Eq. 
However, at low and intermediate velocities substantial changes in the stop- 
ping power are observed as a realistic description of the surface response is 
considered, which leads to the conclusion that a self-consistent description 
of the surface response is necessary if one is to look at the energy loss of 
charged particles moving outside a solid surface. Accurate measurements of 
the energy loss of protons being reflected from a variety of solid surfaces at 
grazing incidence have been reported [THl 113 EO] • A theoretical description 
of these experiments requires that the ion trajectory z(t) be calculated and 
energy losses from the excitation of inner shells be taken into account. Work 
in this direction is now in progress. 
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Finally, we have extended existing first-principles calculations of the ran- 
dom stopping power of valence electrons in Al and Si by allowing transitions 
to larger values of the momentum transfer. We have found that at low veloci- 
ties (where losses from the excitation of inner shells is negligible) the random 
stopping power of Si is considerably smaller than that of Al, though both 
Al and Si have nearly the same valence-electron density. At high velocities, 
band-structure effects become negligible and the random stopping power of 
valence electrons in Al and Si nearly coincide. 

A quantitative comparison of our theory with existing measurements of 
the energy loss of antiprotons jT2] (which unlike protons carry no bound 
states) in a variety of target materials can be achieved by combining our 
first-principles calculations of the Z\ (linear-response) stopping power with 
Z\ corrections in a FEG. Nevertheless, a comparison with experiment still 
requires the inclusion of losses from the inner shells, xc effects, and higher- 
order nonlinear terms. Work in this direction is now in progress. 
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